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■ ABSTRACT 

We present results of a statistical analysis of the SFI catalog of peculiar velocities, a recently completed 
survey of spiral field galaxies with l-band TuUy-Fisher distances. The velocity field statistic utilized is 

■ the velocity correlation function, ipi{r), originally introduced by Gorski et al. (1989). The analysis is 
I performed in redshift space, so as to circumvent potential ambiguities connected with inhomogeneous 

Malmquist bias corrections. The results from the SFI sample are compared with linear-theory pre- 
dictions for a class of cosmological models. We generate a large set of mock samples, extracted from 
N-body simulations, which are used to assess the reliability of our analysis and to estimate the associ- 
ated uncertainties. We assume a class of CDM-like power spectrum models, specified by erg, the r.m.s. 
fluctuation amplitude within a sphere of 8/i^^Mpc radius, and by the shape parameter T. Defining 
?78 = (Tsrio'^j we find that the measured V'i('') implies a degenerate constraint in the r/s-T plane, with 
rjs = 0.3 ± 0.1(r/0.2)°-^, at the 2a level, for the inverse TuUy-Fisher (ITF) calibration presented in this 
paper. We investigate by how much this constraint changes as we account for uncertainties in the anal- 
ysis method and uncertainties in the distance indicator, and consider alternative ITF calibrations. We 
find that both changing the error weighting scheme and selecting galaxies according to different limiting 
line-widths has a negligible effect. On the contrary, the model constraints are quite sensitive to the ITF 
calibration. The other ITF calibrations by Giovanelli et al. (1997) and da Costa et al. (1998) both give, 
for r = 0.2, r]s c:^ 0.6 as the best-fitting value. FERMILAB-Pub-99/133-A 
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1 INTRODUCTION 

The peculiar velocity field of galaxies provides a very 
powerful way of probing mass fluctuations on intermedi- 
ate to large scales ( < 100 Mpc, h being the Hubble 
constant in units of lOOkms"^ Mpc~^), as it is sensitive 
primarily to large scale density fluctuations. Therefore, 
studies of cosmic flows can be \iscd to constrain the ampli- 
tude of the large-scale mass power-spectrum, thus comple- 
menting the information on intermediate scales, between 
those probed by rcdshift surveys and those sampled by 
anisotropics in the cosmic microwave background (CMB) 
as observed by COBE (see the review by Dekel 1994). An- 
other advantage in studying the velocity field is that it is 
measured on scales where linear approximation to gravita- 
tional instability is expected to hold, thus allowing one to 
explore more thoroughly the parameter space of cosmolog- 
ical models. We can parameterize the fluctuation power 
spectrum in terms of the r.m.s. fluctuation within spheres 
of 8h~^Mpc, CTg, and of a shape parameter T. Then, ac- 
cording to linear theory, the typical amplitude of the pecu- 
liar velocity on a given scale is proportional to % f{r, R), 
where t/8 = (Jsf^m^ (following the notation of Chiu, Os- 
triker & Strauss 1998; here is the matter density pa- 
rameter) and f{T,R) is a quantity which depends on the 
power spectrum shape and on the scale R at which the 
velocity field is probed. 

Several statistical characterizations of the peculiar 
velocity fields have been proposed in the last decade, with 
the aim of providing more robust constraints on cosmolog- 
ical scenarios, as newer and larger data sets came to com- 
pletion (e.g., Strauss & Willick 1995, for a review). Among 
such statistical measures, in this paper we will concentrate 
on the velocity correlation function, which has been intro- 
duced for turbulence studies by Monin & Yaglom (1975) 
and borrowed for cosmology by Peebles (1980; cf. also 
Gorski 1988). We will apply this statistic to the SFI sam- 
ple, a recently completed homogeneous all-sky survey of 
Sbc-Sc galaxies with I-band TuUy-Fisher (TF) distances 
(Giovanein et al. 1997a; Haynes et al. 1999a,b, H99a,b). 

A first application of the velocity correlation statis- 
tics to observational data was realized by Gorski et al. 
(1989, G89 hereafter; see also Groth, Juszkiewicz & Os- 
triker 1989), who analyzed the spiral galaxy sample by 
Aaronson, Huchra & Mould (1979) and the elliptical 
galaxy sample by Burstein et al. (1987), finding substan- 
tial discrepancies between the results obtained from these 
two data sets. Tormen et al. (1993, T93) analyzed the 
correlation statistics of the Mark II sample, with results 
favoring 773 ~ 0.7 for scale-invariant CDM models. Kolatt 
& Dekel (1996) estimated the matter power-spectrum im- 
plied by the POTENT reconstruction of the Mark III data 
(Willick et al. 1997) and found % 0.7-0.8. More re- 
cently, maximum-likelihood analyses, estimating the mass 
power-spectrum that gives rise to the observed peculiar 



velocities, have been performed by Zaroubi et al. (1997) 
on the Mark III sample and by Freudling et al. (1999, 
FZ99) on the SFI sample. Both analyses consistently 
find ?78 ~ 0.8 ± 0.2 (90% c.L), quite independent of the 
power-spectrum shape. These results point toward high- 
amplitude fluctuations, thus somewhat at variance with 
results from the r.m.s. cluster peculiar velocity (e.g., Bor- 
gani et al. 1997; Watkins 1997) and with constraints from 
the local cluster abundance (e.g.. Eke, Cole & Frenk 1996; 
Girardi et al. 1998), which indicate lower values. 

Studies of the peculiar velocity can also be combined 

with analyses of all sky rcdshift surveys to investigate the 
relation between the galaxy and underlying mass distribu- 
tions, a key ingredient for understanding galaxy biasing. 
Comparisons between the measured peculiar velocities or 
the recovered densities with those predicted from all-sky 
redshift surveys are commonly used to estimate the pa- 
rameter /3 = fi°„®/6, under the assumption of linear bias- 
ing with a bias factor b. Several estimates of /? have been 
presented in the literature (e.g., da Costa et al. 1998; 
Wilhck & Strauss 1998; Branchini et al. 1999, and refer- 
ences therein) based on comparisons between the velocity 
flelds directly inferred from TF data and recovered from 
galaxy density field in the IRAS 1.2 Jy (Fisher et al. 1995) 
and PSCz survey. Such analyses generally find /? values in 
the range 0.5-0.7. Taking b = (J%,iras/<^9,: these results 
would imply 778 ~ 0.35-0.50 for (Jsjras — 0.7 (Fisher et 
al. 1994). On the other hand, analyses based on the com- 
parison of density fields provide values of /3 as large as 
0.9 (e.g., Sigad et al. 1998). The interpretation of the 
(3 values is further complicated if galaxy biasing is better 
described by a stochastic, nonlinear process (e.g., Dekel & 
Lahav 1999). 

The aim of this paper is to perform a detailed anal- 
ysis of the velocity correlation function for the SFI sam- 
ple and to derive the resulting constraints on large-scale 
structure formation models. The comparison to theoreti- 
cal expectations is based on linear-theory predictions and 
we resort to large-scale N-body simulations to verify the 
reliability of our analysis and to estimate the associated 
errors, contributed by both the cosmic variance and by 
the scatter in the TF relation. 

In our analysis, we choose to use rcdshift space in- 
formation as the indicator of distance for the SFI galaxies, 
so as to avoid the associated Malmquist bias arising from 
the intrinsic scatter of the distance indicator when using 
the inferred distances (cf. Freudling et al. 1995, for a 
discussion on bias corrections in the SFI sample). The 
forward TF relation, obtained by regressing the apparent 
magnitudes over the line-width, in this case, is still suscep- 
tible to selection bias due to the imposed magnitude-limit. 
Using the inverse relation, i.e. fitting the line width as a 
function of the apparent magnitude, avoids this selection 
bias, as long as the sample selection is independent of the 
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line-width (see §6 of Strauss & Willick 1995, and refer- 
ences therein). For this reason, we perform our analysis in 
redshift-space by using peculiar velocities estimated from 
the inverse TuUy-Fisher (ITF) relation. 

The outline of the paper is as follows. In Section 2 we 
provide a basic description of the SFI sample and present 
the ITF calibrations on which our analysis is based. Sec- 
tion 3 contains a brief introduction to the velocity correla- 
tion formalism and presents the results of its application 
to the SFI data. In Section 4 we present the velocity corre- 
lation analysis of our mock samples. In Section 5 we derive 
the resulting constraints on cosmological models and dis- 
cuss the impact of systematic effects in both the sample 
definition and the correlation analysis method. We sum- 
marize our main conclusions in Section 6. 



2 THE SFI SAMPLE 

The TF data defining the sample used here consists 
of two main datasets: a subset of the Mathewson, Ford 
& Buchhorn (1992) survey with about 1200 galaxies with 
I-band photometry and measured rotational velocities, ei- 
ther from radio observations of 21-cm line-widths or opti- 
cal rotation curves; the SFI I-band TF redshift-distance 
survey of about 1300 Sbc-Sc field galaxies. The SFI sam- 
ple consists of galaxies with inclination > 45° north of 
6 < —45° and galactic latitudes \b\ > 10°. The original 
Mathewson et al. (1992) measurements of magnitude and 
rotational velocities were converted into the SFI system 
using about 200 to 300 common galaxies. 

In addition to the field galaxies, roughly 800 galaxies 
covering a broader range of morphological types were ob- 
served in the field of 24 clusters (Giovanelli et al. 1997a, b; 
SCI sample). After careful membership assignment, clus- 
ter galaxies were used to derive a combined TF relation 
corrected for Malmquist bias and bias introduced by in- 
completeness and different morphological mix. In order 
to perform our analysis in redshift-space, we consider the 
inverse TF relation (ITF, hereafter) between the absolute 
magnitude M and the full line width W, 

M = a + b{logW - 2.5) , (1) 

with a = —20.95 and b = —7.94 (here W is expressed 
in units of km s^^ and we assume a Hubble constant of 
100 km Mpc~^). This relation has the same slope as 
that originally provided by Giovanelli et al. (1997b, G97 
hereafter), whose zero-point, a = —21.10, is 0.15 mag- 
nitudes smaller. This difference is due to a new deter- 
mination of the velocity widths and to the removal of 71 
galaxies due to poor photometry, poor line-widths or ob- 
vious misidentification (cf. H99a,b). The la uncertainty 
in the zero-point has been estimated by G97 to be about 
0.05 magnitudes, when combining statistical uncertainties 



in the TF fitting and uncertainties in defining the cluster 
reference frame with a finite number (24) of such objects. 
This uncertainty does not however include possible sys- 
tematics associated with the processing of the raw data 
or with difference between the TF relation of clusters and 
field galaxies, or potential deviations of our local universe 
from a global Hubble flow (e.g. Zehavi et al. 1998, but see 
also Giovanelli et al. 1999). 

We note that careful analysis of the TF relation for 
galaxies in clusters suggests that the scatter depends on 
the line-width. This dependence is modeled by letting the 
error in the estimated distance of the i- th galaxy to 
be Cj = A{Wi)ri, where A(Wi) is the fractional error in 
the distance as estimated from the scatter about the ITF 
relation as a function of the measured line-width of the 
galaxy (G97, cf. also Willick et al. 1997 and Willick & 
Strauss 1998). The resulting errors are estimated to be in 
the range 15-20%. 

Unless otherwise specified and following da Costa et 

al. (1996) and FZ99, we discard those (- 7%) SFI galaxies 
with line-width logVF < 2.25, because of the limited relia- 
bility of the ITF relation at such line-widths. We will also 
show the robustness of the final results against changes in 
the assumed limiting line-width. Furthermore, we restrict 
our analysis to the SFI subsample defined by galaxies ly- 
ing within cz < 6000 km s"'^. With such restrictions, the 
final sample on which we base our analysis contains 974 
galaxies. 

A further alternative calibration of the ITF has been 
presented by da Costa et al. (1998, dC98 hereafter), based 
on a comparison of the velocity field of the SFI sample 
and that implied by the IRAS 1.2 Jy survey. The result- 
ing zero-point and slope of the ITF are a = —21.11 and 
b = —8.55, respectively. In the following, we will use the 
above most recent ITF calibration as the reference one, 
but will show the eSect of taking the previous G97 and 
dC98 calibrations on the final constraints on cosmological 
parameters. 



3 THE VELOCITY CORRELATION 
STATISTICS 

The estimator for the velocity correlations that we 
will use in the following is that introduced by G89 and is 

given by 

Mr) = Yl ' 

where 'dij is the angle between the direction of the i-th 
and the j-th galaxy and the sums are over all the galaxy 
pairs at separation r in rcdshift space. With the above def- 
inition, the ^'i(r) statistics is independent of any assump- 
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tions regarding the velocity field, such as homogeneity and 
isotropy, and has been shown by G89 to be rather robust 
to sampling fluctuations. In eq.(||) Ui is the radial pecu- 
liar velocity of the i-th galaxy and Wi represents a suitable 
weight to be assigned to it. The introduction of the weights 
is a slight modification of the expression for -01 provided 
by G89 (see also T93). Different weighting schemes will be 
applied in the following: (1) uniform weighting, Wi = 1; 
(2) weighting galaxies according to their distance-error. 



l/e^; (3) weighting according to 



where aj is the variance of the local velocity field. 

The quantity a/ can be interpreted as a line-of- 
sight velocity dispersion and has been introduced in order 
to model possible non-linearities, which generates small- 
scale random motions within virialized regions. Such mo- 
tions, which would give rise to an uncorrelated veloc- 
ity component, are expected to be relatively unimportant 
for the SFI field galaxies, whose peculiar velocity should 
not be much affected by virial motions. A further pos- 
sible interpretation of cry is an unrecognized distance- 
independent error, which is not accounted for by the 
ITF scatter calibrated by using members of distant clus- 
ters (e.g.. Kaiser 1988). FZ99 checked for such a term 
by having it as a further degree of freedom to be con- 
strained by a maximum likelihood approach and found 
(7/ = 200 ± 120kms~^. When resorting to the weight- 
ing scheme (3), we will take tx/ = 150 kms"^, although 
our final results are essentially insensitive to its choice. 

As for the scheme (1), its main drawback is that it 
assigns the same weight to all objects, regardless of the 
uncertainty in the velocity errors, which increase with dis- 
tance. Although the methods (2) and (3) overcome this 
limitation, they reduce the effective sampling volume, and 
have been shown by Dekel, Bertschinger & Faber (1990) 
to overestimate the contribution of well sampled regions 
with respect to under-sampled regions in the reconstruc- 
tion of velocity fields. In the following we will mainly base 
our analysis on the uniform-weighting scheme, which is 
the least affected by cosmic scatter (see Section ^ below) . 

As shown by G89, the ensemble average of ^pi{r) is 
given by 

*i(r) = (iPiir)) = ^(r)1'||(r) + [l-^(r)]1'i(r), (3) 

under the assumption of homogeneity and isotropy, where 
V&ll and arc the radial and transverse correlation func- 
tions of the three-dimensional peculiar velocity field. In 
linear theory, they are connected to the power-spectrum 
of density fiuctuations, P{k), according to 



*l|(r) = 



27r2 



dkP{k) 
dkP{k) 



jo{kr) 
ji(fcr) 



Jiikr) 



kr 



(4) 



27r2 J kr ' 

where ji{x) is the i-th order spherical Bessel function and 



The quantity A appearing in cq.(|^) is a moment of 
the selection function of the sample depending on the spa- 
tial distribution of galaxies according to 



A{r) = 

X]|ri-rj|=r^i^j [rirj{cos'dij - 1) -l-r^ cosily] cost^j. 
r'^Y,\r,-Vj\=r cos^ 



(5) 



This quantity provides in a sense the relative contribution 
to ipi (r) from the radial and transverse components of the 
velocity correlation. The definition of eq.(||) is slightly dif- 
ferent from that previously adopted by other authors, by 
including the galaxy weights. 

The advantage of using -01 is that it can directly 
be calculated from the observed radial velocities, without 
the need of any additional assumption. It can then be 
related to theory (eq.||],[E|), taking into account the spe- 
cific sampling through eq.^) . The geometrical factor A{r) 
is plotted in Figure 1 for the three mentioned weighting 
schemes. The net effect of a non-uniform weighting is that 
of increasing A{r) at separations ^ 2000 km s^^. This is 
the consequence of the fact that takes relatively more 
contribution than from large-scale fluctuations (see, 
e.g., Gorski 1988). Therefore, its contribution to ^'i(r) 
is suppressed with the error weighting, which amounts to 
decreasing the effective volume of the sample. 
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Fig. 1. — The geometrical factor A{r) (eq.j^; see text), asso- 
ciated with the SFI sample, for the three alternative weighting 
schemes. 
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The velocity correlation function ipi{r) for the SFI 
sample, with the H99 calibration, computed within bins 
of 500kms~^, is plotted in Figure 2. No errorbars are 
assigned here to ipii^)- We will discuss in the next sec- 
tion how to associate uncertainties to model predictions, 
in order to provide confidence levels in the estimate of cos- 
mological parameters. The upper panel shows the effect of 
adopting different weighting schemes. It is apparent that 
the choice for Wi has a marginal impact on the correlation 
signal. This result might seem somewhat unexpected, in 
view of the different A{r) values for the weighted and un- 
weighted cases. However, these differences appear only 
at rather large separations, 2000 kms""'^ (cf. Figure 1), 
where the value of Tpi for SFI rapidly declines, thus making 
any difference among different weighting schemes hardly 
detectable. By comparing this result with that from the 
real-space analysis of the Mark II sample by Tormen et 
al. (1993), it turns out that the SFI sample produces a 
velocity correlation signal which is at least a factor two 
smaller, although the corresponding scales at which V'i(^) 
approaches zero (~ 3000 km s^^) are similar. 



Hubble-like flow er. This change corresponds to the 2a 
formal statistical uncertainty estimated from the analy- 
sis of the SCI sample of cluster galaxies (G97, H99a,b). A 
global Hubble-like flow represents a coherent velocity field 
which is characterized by a positive correlation (i.e., galax- 
ies moving in the same direction) on intermediate scales, 
5000 km s~^ and by a negative correlation at the largest 
scales, 7000 km s~^, when the two galaxies of a pair are 
placed in the opposite directions of the sample. 

Alternative estimators of the velocity correlation 
statistics have been applied by different authors. Groth 
et al. (1989; cf. also Kaiser 1988) considered the generic 
form for the velocity correlation tensor under the assump- 
tion of homogeneous and isotropic velocity field, '^ij{r) — 
{vi{x)vj{x-r)) = '^±{r)6ij + [^\\ir) -'if±{r)]fifj, where 
6ij is the Kronecker symbol. Then, they obtained ^'j, 
and ^Pii by a x^~niiii™ization procedure to the data. 
G89 compared this method to their 5'i(r) approach and 
showed that they produce comparable results, although 
the former turns out to be noisier at large separations, 
r> 4000kms-i. 
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Fig. 2. — The velocity correlation function, ■01 (r) (in units 
of 10*kms~^), for the SFI sample. The upper panel shows the 
effect of different galaxy weights, while the lower panel shows 
the effect of changing by 0.1 magnitudes the zero-point of the 
ITF relation, representing the 2a uncertainty in its calibration 
(cf. G97, H99a,b). 



More recently, Ferreira et al. (1999) proposed a 
new method to estimate the main galaxy pairwise velocity, 
vi2 — {v{xi) — v(x2))- This method, which has been so 
far tested on N-body mock samples and is in the process 
of being applied to real data sets, provides essentially con- 
straints on CT^rJ^^. Therefore, its combination with linear- 
theory constraints on us^^m^ could in principle break the 
degeneracy between erg and fi™. Of course, careful in- 
vestigations are required in order to understand whether 
available data are of sufficient quality and their systemat- 
ics and biases are enough under control to allow a reliable 
estimate of cg and separately. 



4 ANALYSIS OF THE MOCK SAMPLES 



In order to explore extensively the model parame- 
ter space, we resort in the following to linear theory as the 
means to compare model predictions and SFI results. Two 
important issues need to be addressed: (a) the reliability 
of our analysis and specifically the use of linear theory 
to predict the statistics of the velocity field, and (6) the 
estimate of the cosmic scatter and the observational un- 
certainties associated with the SFI sampling, in order to 
establish the confidence level for model exclusion. For this 
purpose we use large N-body simulations from which we 
extract sets of mock samples which mimic the sampling 
and selection effects of the SFI sample. 



The lower panel of Fig. 2 shows the effect of chang- 
ing the zero-point of the ITF relation (eq.j^]) by 0.1 
magnitudes either way, which corresponds to a change 
of e ^ 2.5% in the distances or an additional global 
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4-1 Generating the mock samples 

The parent N-body simulations from which we ex- 
tract mock samples have been run by using the publicly 
available adaptive P^M code by Couchman (1991). We 
have run two simulations corresponding to two different 
cosmological scenarios. The first model is a flat low- 
density one with ftm = 0.4 (AO. 4). The transfer function 
used is that of Bardeen et al. (1986) [see eq. 

(0) below], 

with the shape parameter T set to 0.22 and as = 0.87. 
The second model is an Einstein-de Sitter (EdS) universe, 
with r = 0.43, and erg = 1.2. With the above parameters, 
both models are consistent with the 4-year COBE normal- 
ization (e.g., Bunn & White 1997), while the EdS model 
fails to match the abundance of local galaxy clusters (e.g.. 
Eke et al. 1996; Girardi et al. 1998) and the shape of 
the galaxy power-spectrum (e.g.. Peacock & Dodds 1994; 
Liddle et al. 1996). 

Each simulation follows 128^ particles within a box 
of 250/i~^Mpc on a side. The adopted Plummer softening 
scale, ~ 100 kpc, is more than adequate to describe 
the large-scale velocity field (see Borgani et al. 1999, for 
a more detailed description of the simulations). Velocity 
fields on scales of a fewx 10/i^^Mpc, which are of interest 
in this paper, receive a small but non-negligible contri- 
bution from wavelengths larger than the adopted box size. 
Furthermore, the volume of a single simulation can accom- 
modate only a rather small number of non-overlapping SFI 
mock samples (each extending out to cz — 6000 km s^^), 
so as to not allow a reliable determination of cosmic vari- 
ance. 

In order to extend the dynamic range of our simula- 
tions to larger scales, we resorted to the method proposed 
by Tormen & Bertschinger (1996) of adding longer waves 
to N-body outputs. This method, which allows to gen- 
erate non-periodic replicas of a parent box, is based on 
the Zel'dovich approximation (Zel'dovich 1970) for com- 
puting the contribution to particle displacements and ve- 
locities from waves longer than the original box size. Cole 
(1997) showed that this procedure is adequate to extend to 
larger scales the description of peculiar velocities. In our 
analysis, we replicate the original box three times along 
each spatial direction, which leads to a total of 27 replica 
and a final box of size L = 750/i~^Mpc, containing about 
5.7 X 10^ particles. 

As a first step for mock sample extraction, we di- 
vide the large box into 6'^ smaller boxes of 125/i~^Mpc 
on a side. At the center of each of them we place an 
observer. After randomly choosing the orientation of the 
"galactic" coordinate system, we select among the simu- 
lation particles those which are closest to the position of 
real galaxies in the SFI sample. In this way, we gener- 
ate mock samples with the same spatial distribution and 
number of galaxies as in the real SFI sample. The "true" 



radial velocities in the mock samples are perturbed ac- 
cording to the associated observational errors of the real 
catalog and according to the assumed random velocity dis- 
persion af (under the assumption that both contributions 
are independent Gaussian variables). For each simulation, 
we generate two sets of mock samples, based on assuming 
both (7/ = and 150 km s""'^. Since the final results turn 
out to be essentially indistinguishable, we will present for 
the mock sample analysis only results based on assuming 
a vanishing a f . 

We note that other authors (e.g., G89; Strauss, Gen 
& Ostriker 1993; T93) followed more sophisticated proce- 
dures to search for "observers" within simulations. Such 
procedures involve selecting observers so that local prop- 
erties of the density and velocity field resemble those ob- 
served for the Local Group of galaxies. However, T93 
showed that applying such constraints on the observer se- 
lection does not significantly alter the velocity correlation 
statistics for realistic power spectra. Furthermore, the aim 
of our analysis is to estimate how often the SFI correlation 
statistics can be observed in a given cosmology assuming 
the variety of observers' characteristics to be included into 
the cosmic variance which is appropriate for that model. 
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Fig. 3. — The comparison between linear-theory predictions 
(dashed curves) and results from the analysis of mock sam- 
ples for the velocity correlation function i/'i(r) (in units of 
10* kms~^). Mock samples are extracted from an N-body sim- 
ulation of the AO. 4 model. Upper and lower panels refer to 
uniform weighting and distance-error weighting, respectively. 
Error bars are the la scatter among the set of 216 mock sam- 
ples. 
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4-2 Testing the analysis method 

Since the mock samples have been generated by re- 
producing the positions of real galaxies, their correspond- 
ing A{r) is the same as for the real SFI sample. For each 
cosmological model we compute in linear theory the ex- 
pected tpi (eq.||]) and compare it to the distribution of 
values obtained from the mock samples using eq.(||). We 
plot in Figure 3 the results of this comparison for the AO. 4 
case, for both uniform (upper panel) and distance-error 
(lower panel) weightings. Filled circles represents ipi (r) as 
estimated by averaging over the set of Nmock = 216 model 
samples and the errorbars are the la scatter, arising from 
both cosmic variance and observational uncertainties. As a 
basic result, it turns out that, for both weighting schemes, 
linear theory is always adequate to describe the expected 
velocity correlation function for samples having the same 
selection effects as the SFI, once they are accounted for 
by the A{r) quantity. Any residual discrepancy on small 
(^ 1500 km s^^) scales, which are probably due to sam- 
pling effects or to residual non-linearities, are well within 
the la scatter. Furthermore, we also remind that since 
the SFI sample only contains field spirals, we expect their 
dynamics to be even closer to linear theory that that of 
the N-body particles belonging to the mock samples, that 
we did not attempt to select so as to avoid high-density 
regions. 

We checked the relative contribution to the errors 
from the cosmic scatter and from the uncertainties in the 
peculiar velocity measurements, using a set of mock sam- 
ples where peculiar velocities are not perturbed according 
to ITF distance errors, so that only the effect of the cos- 
mic scatter is present. It turns out that the cosmic scat- 
ter is clearly dominant at r < 3500 km s~^, with the TF 
scatter contributing < 20% and becoming relevant only at 
larger scales. The distance-error weighting scheme gener- 
ates a larger scatter, as a consequence of the fact that this 
method amounts to reducing the effective volume where 
ipi^r) is computed. For this reason, in the following we 
will take the uniform weighting as the reference analysis 
method to constrain model parameters. 

4-3 Estimating ijji uncertainties 

Having demonstrated that linear theory provides re- 
liable predictions for tpi, the next information that one 
needs is the uncertainty to be associated to such predic- 
tions. In order to do so, we estimate from the set of mock 
samples the elements of the covariance matrix, C*-' , which 
are defined as 

c'' = — E i^h - ^0 [Ki - i^i) ■ (6) 

mock i — \ 

Here ■0! ; is the value of the velocity correlation function 
at the z-th separation bin for the ^-th mock sample, while 



is its average value estimated over the Nmock samples. 

Figure 4 shows the comparison between results from 
the AO. 4 and EdS models, by plotting the quantities 
C^^ /^{tpl. According to its definition, this quantity de- 
scribes the relative covariance of the ipi values at different 
separations. The top panels show the results for the di- 
agonal (variance) terms, while the other panels show the 
off-diagonal terms, illustrating different rows in the co- 
variance matrix. The first thing to note is the large cross- 
correlation between the results of the different bins, which 
are comparable to the variances, and therefore cannot be 
ignored when using the ipi statistic to constrain cosmolog- 
ical models. 

In addition, it is apparent from the figure that, apart 
from small differences due to statistical fiuctuations, the 
two models have the same amount of relative covariance. 
This is not unexpected since, to a first approximation, the 
long-wave perturbations which generate the cosmic scat- 
ter, are also responsible for the ■0i signal, so as to make 
the relative scatter fairly constant. Noticeable differences 
occur only at relatively large separations, > 3500 km s~^, 
where the observational uncertainties become more domi- 
nant, thus increasing the total scatter and suppressing the 
discriminative power of ipi{r)- For this reason, in the fol- 
lowing we will compare linear-theory predictions and SFI 
results only for r < 3500 km s~^, where the relative un- 
certainties are essentially the same for the two considered 
models. We note that, since AO. 4 and EdS have rather 
different values for both ryg and for the power-spectrum 
shape, we can quite confidently conclude that the relative 
scatter for ipiir) is model-independent, at least for the 
range of models and scales of interest, while its absolute 
value is not. 

In the top right panel of Fig. 4 we compare the 
diagonal terms for the AO. 4 mock samples for ipi com- 
puted according to uniform and distance-error weighting 
schemes. It is apparent that the distance-error weighting 
is associated with larger error bars, as was already shown 
in Figure 3. 

Based on these results we, therefore, conclude that: 
(a) the errors of individual "01 bins are significantly cor- 
related; (b) a general recipe can be devised for the ipi 
uncertainties, whose relative amount is fairly independent 
of the cosmological model; and (c) that the size of such 
errors is smaller when ipi is estimated according to the 
uniform- weighting scheme. 



5 CONSTRAINING COSMOLOGICAL MODELS 

Based on the results obtained so far, we will now use 
eqs.(^) and (^ as a model prediction for ipi. As for the 
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Fig. 4. — The elements of the relative covariance matrix, /'ip\%l)\, for SFI mock samples extracted from EdS and AO. 4 simulations. 
The top panels shows the diagonal (variance) terms, with the top right panel comparing the variance for unweighted and error- 
weighted estimates of ipiir). The other panels are for the off-diagonal terms and show different rows in the covariance matrix (see 
text). 



model power spectrum, we express it as P{k) ~ AkT'^{k) 
where we assume a Harrison-Zel'dovich shape on large 
scales. The transfer function, T{k), is taken to be 



nq) 



ln(l + 2.Mq) 
2.34g 



{bA&qf + {&.nqf 



-1/4 



(7) 



where q — k/Th and T is the so-called shape parameter. 
For r ~ flmh, eq.(0) provides the transfer function for 
CDM models with a negligible baryon fraction (Bardcen 
et al. 1986). More generally, it can be seen as a phe- 
nomenological expression, with F a parameter to be fixed 
by observational constraints. As for the amplitude of the 
power spectrum, it is customary to express it in terms 
of as- Following eqs.(Q), the velocity correlation hmction 
ipi (r) is then entirely specified in linear theory by the two 



parameters F and 778. 

Despite the errorbars being so large such that the 
i/'i detection is only marginally different from zero in each 
individual bin (cf. Figs. 2 and 3), its determination at dif- 
ferent scales does allow to place significant constraints on 
the rjs-T plane. In order to provide constraints on these 
parameters, we compute the weighted between the SFI 
correlation function, ipi^^ : that from model predic- 
tions, 



X 



-1 



[i^r\n)-i;r'ir^)] . (8) 

dements of the inverse 
ance matrix, as calibrated from the mock samples, and 



Here, C^j^ are the elements of the inverse of the covari- 
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the sums are over the radial bins of 500 km s"-'^ width, for 
separations r < 3500 km s~^. The probability for model 
rejection is estimated by assuming a statistic, from the 
value of Ax^ = - Xmim where Xmin is the absolute 
minimum value. 




Fig. 5. — The la and 2a contours in the 77g-r plane from the 
analysis of the velocity correlation function, tpi (r) , for different 
calibrations of the inverse TuUy-Fisher relation. The horizontal 
shaded area corresponds to the 90% confidence level constraints 
on ris from the analysis of the Giovanelli et al. (1997a,b) r.m.s. 
cluster peculiar velocities (Borgani et al. 1997). The vertical 
shaded area is the 95% confidence level constraint on the shape 
parameter from the power-spectrum of APM galaxies (Liddle 
et al. 1996). 

In Figure 5 we plot the iso-Ax^ contours for the 
three ITF calibrations of the SFI sample that were dis- 
cussed in Section 2. Internal and external contours corre- 
spond to Ax^ = 2.30 and 6.17, respectively, thus providing 
the 1(7 and 2a confidence levels for two significant param- 
eters. The corresponding minimvim values of the x^ per 



degree of freedom are 1.67, 0.80, and 0.78, for the H99, 
G97 and dC98 calibrations, respectively. In all cases, the 
best-fitting model seems to provide an acceptable fit. This 
value for the H99 calibration is somewhat large, however it 
corresponds to only ^ la deviation for a x^ statistic with 
five degrees of freedom. The fact that such x^ values are 
around unity indicates that our error model is realistic. 

The vertical shaded areas represent the 95% confi- 
dence level interval on the shape parameter, as derived by 
Liddle et al. (1996) from the power-spectrum of APM 
galaxies. The horizontal shaded areas represent the 90% 
confidence level on rjs derived by Borgani et al. (1997) 
from an analysis of the r.m.s. peculiar velocity of SCI 
clusters (Giovanelli et al. 1997a). All these constraints 
intersect our 2a confidence regions. 

For the II99 and G97 ITF calibrations, the con- 
straints in the 778^r plane can be cast in the form 

/ F V 

Vs = Vs,o X ( — 1 , (9) 

with r]s,o = 0.30+°;j7 and ?78,o = 0.58^0^1 for the two 
above calibrations, respectively (errorbars correspond to 
2a c.l.). The asymmetry in the errors is due to the fact 
that, as r]s is increased from its best-fitting values, larger 
absolute errors are assigned to ipi , since the relative scatter 
is taken to be constant (cf. §4.3). Thus, larger values of 
rjg tend to be excluded at a lower significance than smaller 
values. As for the dC98 calibration, the corresponding 
constraints show a somewhat steeper F-dependence of rjg 
with values of F^ 0.35 ruled out at about 2a c.l. It is in- 
teresting to note that, for F ~ 0.2, this result agrees with 
(3 = 0.6 ± 0.1, as found by da Costa et al. (1998), for an 
almost unbiased IRAS galaxy distribution. 

We show in Figure 6 the variation of Ax^ around its 
minimum as a function of yyg, in order to show the effect 
of changing other assumptions underlying our analysis. In 
all the panels, the solid curve refers to constraints from the 
H99 ITF calibration, for a fixed shape-parameter F = 0.2 
and logVF > 2.25 for the line-width of SFI galaxies. 

As demonstrated already (Fig. 2, top panel), our 
results are insensitive to the choice of galaxy weighting, 
and we adopt here throughout the uniform weighting. As 
is illustrated here [panel (6)], changing the limiting line- 
width of the sample also has a negligible effect on our re- 
sults which are virtually unchanged as we increase it from 
2.25 to 2.40. We find as well that our constraints do not 
depend on the specific choice of binning used in the com- 
putation of ipi{r). The effect of the zero-point uncertainty 
is shown in panel (a). As was illustrated also in the lower 
panel of Figure 2, the results are quite sensitive to such 
changes, and a negative shift of the ITF zero-point by 0.1 
mag leads to a sizeable increase of r/g from ~ 0.3 to ~ 0.55. 
For higher values of F, this change would similarly corre- 
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Fig. 6. — Variation of Ax^ around its minimum value as a function of tjs. In both panels the solid curve corresponds to F = 0.2, 
uniform weighting in the estimate of ^l^iir), ITF calibration by H99a,b, with the best-fitting value of the zero-point, and log W > 2.25 
for the galaxy line-width. Panel (a): effect of changing the zero-point of the ITF relation; short- and long-dashed lines are for 
shifting it by 0.1 magnitudes upwards and downwards, respectively. Panel (b): effect of increasing the limiting line-width; short- 
and long-dashed lines are for log W > 2.3 and log W > 2.4, respectively. 



spond to higher values of 778, e.g for T = 0.4, rjs would 
increase from ~ 0.4 to ~ 0.8, and its effect is generally 
comparable to that of varying the ITF calibration. 

Despite the fact that the constraints on cosmologi- 
cal parameters drawn from the V'l statistics are quite sen- 
sitive to the details of the ITF calibrations, some con- 
clusions can still be drawn. First, the constraints on the 
velocity power-spectrum normalization, rjs , depend on the 
P{k) shape, as a consequence of the fact that we are prob- 
ing velocity fields on scales larger that the 8/i~^Mpc nor- 
malization scale. Second, assuming F ~ 0.2, as indicated 
by galaxy clustering data, implies power-spectrum ampli- 
tudes which can be different by up to a factor two, but are 
still generally consistent with independent observational 
constraints. For instance, the local abundance of galaxy 
clusters to a first approximation also provides a constraint 
on = 0.5-0.6 (e.g. Eke, Cole & Frenk 1996; Girardi et 
al. 1998, and references therein). 

Our results for ?7g can also be compared with those 
obtained by Zaroubi et al. (1997) and FZ99, who esti- 
mated the mass power spectrum by a maximum likeli- 
hood (ML) analysis of the peculiar velocities of the Mark 
III and the SFI samples, respectively. These estimates 
are then translated to constraints on 773 by integrating 
over the corresponding spectra. Both works consistently 
found 778 — 0.8 ± 0.2 at 90% c.l. and a preferred value of 
r ~ 0.4 ± 0.2. As the application of the ML analysis for 
the SFI sample has been performed using the G97 calibra- 
tion, it is most suitable to compare the FZ99 results with 



those reported in the central panel of Fig. 5. It turns out 
that the confidence regions coming from the ML and 7/;i 
analyses do overlap over a significant portion of the T^g-F 
plane. For F = 0.4, the tpi analysis gives t]s = 0.85lg :[q. 
The main difference being the dependence of the 778 con- 
strains on F in the ipi analysis, such that for lower values 
of F ~ 0.2 the preferred rjs values are somewhat smaller 
than those obtained in the ML analysis. 

One should also bear in mind the different sensitiv- 
ities of these two analyses. As demonstrated in Figs. 5 
and 6, the ijji analysis is sensitive to the ITF calibration, 
while it is robust to changing the limiting line-width. On 
the other hand, the ML analysis is remarkably robust to 
changes in TF calibration (e.g. Fig. 8 in FZ99), while 
it is more sensitive to the pruning of SFI galaxies at dif- 
ferent line- widths. For these reasons, these two methods 
should be regarded as complementary and both worth to 
be appHed to a given data set. 



6 CONCLUSION 

In this paper we have presented an analysis of the 
velocity correlation function, ipi{r), for the SFI sample 
of Sbc-Sc galaxy peculiar velocities based on the infrared 
TF distance indicator calibrated using a sample of clus- 
ter galaxies (Giovanelli et al. 1997a, b; Haynes et al. 
1999a,b). In order to minimize uncertainties related to 
Malmquist bias corrections, we performed the analysis us- 
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ing the redshift-space positions of galaxies and the ITF 
distance indicator. Three different ITF caUbrations for 
the SFI sample have been examined in our analysis: one 
based on an updated version of the SFI sample presented 
by Haynes et al. (1999a, b, II99), that presented by Gio- 
vanelli et al. (1997b, G97) and that obtained by da Costa 
et al. (1998, dC98). 

The final goal of our analysis is to place constraints 
on the amplitude and the shape of the fluctuation power- 
spectrum, by comparing ipi (r) from SFI and from linear- 
theory predictions of cosmological models. For this pur- 
pose, we needed to verify the reliability of linear -theory to 
predict il^i{r) for a sample having the same galaxy posi- 
tions and observational uncertainties as the SFI one, and 
to estimate the associated uncertainties due to cosmic scat- 
ter and observational uncertainties. These two goals have 
been achieved by comparing linear-theory predictions to 
results from the analysis of a large set of mock SFI sam- 
ples, extracted from N-body simulations. 

We have found that linear-theory provides a rather 
accurate description of the ipi (r) estimated from the mock 
samples, over the whole scale range considered (r < 
5000 km s~^; cf. Figure 3). This confirms that both sparse 
sampling effects and residual non-linearities have a minor 
impact on our analysis. We have also shown that the rel- 
ative covariance in V'i(^) among the set of mock samples 
is roughly independent of the cosmological models, thus 
allowing for a simple treatment of the associated errors. 

In general, we find that our analysis constrains a 
degenerate ridge in the rj^-T plane. For the H99 and 
G97 ITF calibrations, we find r]s = r]s,o(T/0.2f-^ , with 
Vs,o = 0.30_o;o7 and ?78,o = 0.58_g;;^2 ^'^^ the two above 
calibrations, respectively, at the 2a level (cf. Figure 5). 
The dC98 exhibits a stronger tendency for lower values 
of the shape parameter, constraining 0.35 at the 2a 



level, and is consistent with the G97 calibration in that 
range. These constraints are robust to variations of the 
galaxy weighting scheme (cf. Figure 2) and to changes 
in the choice of the limiting galaxy line-width (cf. Fig- 
ure 6), but are, clearly, very sensitive to uncertainties in 
the calibration details, such as the zero-point of the TF 
relation. 

In any case, the results presented here indicate that 
the large-scale velocity field can be brought into agree- 
ment with the low fluctuation amplitude implied at ~ 
10/i^^Mpc scale by the abundance of galaxy clusters (e.g. 
Eke et al. 1996, Girardi et al. 1998), for power-spectrum 
shapes which are consistent with large-scale clustering 
data (e.g. Liddlc et al. 1996), while higher amplitudes are 
allowed for larger values of the shape parameter. Our con- 
straints on the ris-T plane for the ITF G97 calibration and 
those from the maximum-likelihood (ML) analysis for the 
G97 direct TF relation by Freudhng et al. (1999, FZ99) 
are quite consistent for F^ 0.3. Since the ML and the 
tpi methods are sensitive to different degrees to different 
aspects of the analysis (i.e., TF calibration and limiting 
line-width), they should be regarded as complementary 
approaches for extracting cosmological constraints from 
large-scale cosmic flows. 
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